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We analyze the constraints that causality imposes on some of the particle detector models em¬ 
ployed in quantum field theory in general and, in particular, on those used in quantum optics (or 
superconducting circuits) to model atoms interacting with light. Namely, we show that disallowing 
faster-than-light communication can impose severe constraints on the applicability of particle de¬ 
tector models in three different common scenarios: 1) when the detectors are spatially smeared, 2) 
when a UV cutoff is introduced in the theory and 3) under one of the most typical approximations 
made in quantum optics: the rotating-wave approximation. We identify the scenarios in which the 
models’ causal behaviour can and cannot be cured. 


I. INTRODUCTION 

Particle detector models may be thought of as control¬ 
lable quantum systems that couple locally in space and 
time to quantum fields. They provide a way to extract lo¬ 
calized spatio-temporal information from the fields with¬ 
out having to implement projective measurements of lo¬ 
calized field observables mm- This simplifies the task 
of extracting localized information about the field, and 
avoids all the possible complications that may appear 
with the use of projective measurements altogether [^. 
Particle detector models in quantum field theory were pi¬ 
oneered by Unruh and DeWitt ms, and can be found in 
the literature in several slightly different (but fundamen¬ 
tally similar) formats, e.g., a field in a box [i], a two-level 
system [5] or a harmonic oscillator (see, e.g., [SHE]). 

Particle detector models have been successfully em¬ 
ployed in a plethora of contexts in fundamental quantum 
filed theory [a HO]. Perhaps one of the best-known ones 
is the operational formulation of the Hawking and Unruh 
effects (see, e.g., min]). Besides their many uses in fun¬ 
damental quantum field theory, particle detector models 
are ubiquitous as models for experimental setups in quan¬ 
tum optics [12] and in superconducting circuits [T3|. For 
example, an alkali atom as a first quantized system, can 
serve as such a detector for the second quantized electro¬ 
magnetic field. In fact, the common light-matter interac¬ 
tion models, such as for instance the Jaynes-Cummings 
model and its variants, are almost identical in nature 
to the Unruh-DeWitt (UDW) model [12]. Indeed, the 
Unruh-DeWitt detector-field interaction has been proven 
to be a good model of the light-matter interaction in 
quantum optics for processes not involving exchange of 
orbital angular momentum [B US- 

More recently, UDW detectors have been extensively 
used in studies of relativistic quantum information. For 
example, relativistic quantum computing [Bdl] , quan¬ 
tum communication via field quanta |I8II20] . cosmology 
[SHIS, and, more generally, in studies on a host of ef¬ 
fects related to the presence of spacelike entanglement 
in the vacuum state of quantum fields both from funda¬ 
mental [21H?T] and applied [321133] perspectives. Inter¬ 


estingly, in these studies it is shown that it is possible to 
harvest correlations from the field vacuum to spacelike 
separated detectors, which gives an operational proof of 
the spacelike entanglement present in the quantum vac¬ 
uum [531135] . 

Because UDW detector models have been proven so 
valuable in relativistic quantum information, the ques¬ 
tion arises as to what extent these models (which in¬ 
volve non-relativistic systems coupled to fully relativis¬ 
tic quantum fields) behave in a causal way. This is of 
special importance when studying phenomena for which 
the causal behaviour of the model is paramount, such as 
the aforementioned spacelike vacuum entanglement har¬ 
vesting or, more generally, any quantum communication 
scenario where relativistic effects become important. 

It is known that naive projective measurements in 
quantum fields are not compatible with causality and can 
suffer from superluminal signalling (see for instance 13]). 
Do UDW-like detector models have similar problems? 
In particular, in |36] this same question is posed about 
the causal behaviour of the Unruh-Dewitt model. The 
question is well aimed: consider, for example, the usual 
oversimplified single-mode Jaynes-Cummings model m- 
This model allows for faster-than-light signalling and is 
indeed unable to model the light matter interaction in 
relativistic regimes. The question whether the full UDW 
model respects causality was addressed in [T5], where it 
was shown that communication using pointlike Unruh- 
DeWitt detectors as emitters and receivers of field quanta 
is causal to any order in perturbation theory, as long as 
the detectors remain pointlike as originally proposed. 

However, as we will discuss, it is all too common to con¬ 
sider modifications of the original model that sometimes 
jeopardize its causal behaviour, both in theoretical stud¬ 
ies on field theory and in the applications of these mod¬ 
els to quantum optics. Studying the causal behaviour 
of these modified models has a double interest: on the 
one hand, from a theoretician’s point of view, in the sce¬ 
narios where the model presents causality problems it is 
interesting to know if and when the causal behaviour can 
be restored in an approximate way in some regimes. On 
the other hand, as an experimentalist using these models 
to predict phenomenology, it is fundamental to know un- 
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der what circumstances the models can be used to make 
meaningful predictions and, conversely, when the models 
dramatically fail due to violations of a first principle of 
relativistic theory. 

The first modification of the original model whose 
causal behaviour we are going to study is the general¬ 
ization to non-pointlike detectors. It is indeed very com¬ 
mon to consider that the detectors have a spatial smear¬ 
ing [371 EH]. The smearing of particle detectors may re¬ 
spond to the need to regularize divergences of the point¬ 
like model EH]- Smearings can also be introduced, for 
example in quantum optics, to improve on the accuracy 
of the light-matter interaction models assuming that the 
atoms are not pointlike objects, but are instead local¬ 
ized in the full extension of their atomic wavefunctions 
[Ulli]. Furthermore, as discussed in [53] oftentimes spa¬ 
tial smearings are introduced implicitly, hidden in some 
form of soft UV regularization of the model (see, e.g., 

m)- 

As was already noted in [53], and as we will study in 
detail in this paper, spatially compactly supported de¬ 
tectors are safe in terms of their causal behaviour. How¬ 
ever, the use of compactly supported spatial smearing 
is not so common in the literature. Instead, the most 
common smeared UDW models used in the literature 
(both in QFT and quantum optics) assign non-compact 
—but very strongly supported in a finite region— smear¬ 
ing functions, such as Gaussian (See, e.g. [21 EO]) or 
Lorentzian (See, e.g, [57lI59] l smearings. Non-compact 
smearings are employed mainly for two reasons: 1) the 
simplification of analytic computations and the regular¬ 
ization of divergences and, 2) In the context of quan¬ 
tum optics, the introduction of smearing is part of the 
refinement of the model: the atomic wavefunctions of 
valence electrons characterizing the effective size of the 
atom are exponentially suppressed in the radial distance 
to the centre of mass of the atom m- 

A relevant question to ask in these scenarios is to what 
extent the use of fast-decaying smearing functions (as op¬ 
posed to strictly compactly supported ones) effectively 
renders the predictions of the model acausal. Rigor¬ 
ously speaking, any non-compactly supported smearing 
allows for instantaneous signalling. For example, if we 
consider communication between two non-compact de¬ 
tectors A and B, the influence of A in B would instan¬ 
taneously be felt by B. Nevertheless, if in spite of being 
non-compact, the smearing of the two detectors is only 
strongly supported around a point of space with a charac¬ 
teristic length a (as is the case of Gaussian or Lorentzian 
smearing), causality may be recovered in an approximate 
sense when we consider, for instance, communication be¬ 
tween two detectors whose spatial separation is much 
larger than a. In this paper we will study quantitatively 
how the leading-order causality violations in the commu¬ 
nication between two detectors are suppressed faster than 
exponentially as the characteristic size-scale of the non¬ 
compact smearing is reduced. Hence, this leads us to con¬ 
clude that the causality of the model can indeed be ap¬ 


proximately recovered with arbitrary precision when us¬ 
ing strongly-supported but non-compact smearings such 
as Gaussian. 

Another very typical modification of the UDW model, 
fairly common in quantum optics, is the introduction of 
hard UV cutoffs in the detector’s proper reference frame. 
This is usually justified by the fact that the atomic re¬ 
sponse to electromagnetic radiation is a function of fre¬ 
quency and it does eventually become negligible in the 
far UV limit. If we were to consider such as an effec¬ 
tive model, we would have to keep in mind that when 
a UV cutoff is introduced the causality of the model 
is again compromised. How causality is impacted by 
this sort of hard UV cutoffs has been studied in cavity 
settings both perturbatively m and non-perturbatively 
[7]. These studies shown that causality violating vanish 
decrease polynomially fast as the UV cutoff is relaxed. 
However, these studies had two limitations: they were 
limited to 1-|-1-dimensional setups and only to a discrete 
set of field modes corresponding to quantum fields in pe¬ 
riodic and Dirichlet cavities. The fact that we consider 
a continuum of modes (field in free space) as opposed 
to a discrete number of them does introduce differences 
in the causal behaviour of the UDW when UV cutoffs 
are introduced. Much more important, as we will see in 
this work, the dimensionality of spacetime is critical for 
the causal behaviour of the communication using UDW 
detectors in the presence of UV cutoffs. For example, 
we will discuss that in a 3-1-1-dimensional scenario the 
causality loss due to the existence of these cutoffs cannot 
be made arbitrarily small just increasing the cutoff with¬ 
out further regularization, in stark contrast with what 
was obtained for 1-|-1-dimensional cavities in mm- 

Finally, we will analyze the impact of causality on 
the validity of the rotating-wave approximation (RWA). 
The RWA [12] is overwhelmingly present in the litera¬ 
ture on quantum optics. The approximation consists of 
the simplification of the light-matter interaction Hamil¬ 
tonian removing the terms that oscillate rapidly. The 
RWA is discussed to be valid for long interaction times 
(much larger than the detector’s characteristic Heisen¬ 
berg time). However, this approximation may result ex¬ 
tremely harmful for the causality of the model. This issue 
has already drawn the attention of several researchers in 
early [inmi] and more recent [21 [53] works. We will an¬ 
alyze the faster-than-light signalling that this approxima¬ 
tion enables in the context of quantum communication of 
particle detectors through the exchange of field quanta. 
We will show that even for long interaction times, the 
acausal signature of this extremely common approxima¬ 
tion is only very slowly erased as the limit of longer times 
is taken, rendering this approximated model unfit to de¬ 
scribe any setting where the causality of the theoretical 
model is important. 

With this perturbative study of the causal behaviour 
of the UDW model and the light-matter interaction, we 
will quantitatively characterize if (and in what regimes) 
the approximations and considerations described above 
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can still be used mantaining some degree of approximate 
causal behaviour. This is important because current ex¬ 
perimental techniques in quantum optics and supercon¬ 
ducting circuits are nearing the point where intrinsically 
relativistic phenomena can be accessed in experiments 
(seee.g. [44]). Finally, we note that we do a general study 
in one, two and three spatial dimensions with a double 
purpose: 1) There exist experimental setups where quan¬ 
tum fields live in an effectively reduced dimensionality 
such as superconducting microwave guides coupling to 
artificial atoms nmsiiisj, and 2) in doing so, we will 
be able to characterize the relevant role that the number 
of spatial dimensions has in the causal behaviour of the 
light-matter interaction models, and learn some interest¬ 
ing aspects of the regular UV behaviour of the signalling 
between two particle detectors. 

II. SIGNALLING OF TWO PARTIGLE 
DETECTORS THROUGH A QUANTUM FIELD 

We will study the Unruh-DeWitt particle detector 
(from now on referred to as the ‘detector’) which will 
be considered in this work as a spatially localized two- 
level quantum system coupled to a scalar field. The spa¬ 
tial smearing of the detector will be given by the real¬ 
valued smearing function F{x) which we choose, for con¬ 
venience and w.l.g., to be localized and centred around 
a; = 0. This detector couples to a scalar field locally 
along its trajectory in spacetime. In general, we can con¬ 
sider that the detector’s centre of mass moves in a trajec¬ 
tory parametrized in terms of its proper time r. That is, 
X = Xq{t). In the particular case of a stationary detector 
whose centre of mass is placed at a; = a;o, the detector is 
comoving with the usual Minkowski frame (a:, t) in which 
we carry out the field quantization. This means that we 
can take t = t. The Unruh-DeWitt interaction Hamilto¬ 
nian (in the interaction picture) for this case is given by 

m 

Hi = X J d'^x F[x — Xo]x{t)m{t)^[x,t] (I) 

where x(t) is the switching function controlling the 
coupling-decoupling speed and the duration of the 
detector-field interaction, m{t) = (d+e’^* -I- is 

the detector’s monopole moment (being Q the energy gap 
between the detector’s two energy levels), ^[x,t] is the 
quantum scalar field, and n is the number of spatial di¬ 
mensions. The difference of 0 with the pointlike model 
is that the field operator is evaluated along the worldline 
of the detector, as usual, but summed over the whole 
spatial extension of the detector. 

Although this is out of the scope of the present pa¬ 
per, it is worth mentioning that if we were to generalize 
the analysis of these localized detector models to general 
non-inertial trajectories, we would have to face the well- 
known problem of accelerating rigid bodies. Roughly 
speaking, the proper distance between two points of a 


solid accelerating with the same proper acceleration in¬ 
creases in time, eventually destroying the solid when the 
internal cohesion forces that support it are overrun by 
“stress forces” m- A reasonable hypothesis made for a 
physical detector in those cases is that, until this hap¬ 
pens, the detector has to keep internal coherence to a 
good approximation. In other words, the internal forces 
that keep the detector together will prevent it from be¬ 
ing further smeared due to relativistic effects, up to some 
threshold acceleration. This in turn implies that, effec¬ 
tively, every point of the detector will accelerate with a 
different acceleration in order to keep up with the rest of 
its points and maintain the shape of the spatial profile of 
the detector in the centre-of-mass reference frame. The 
natural framework to treat such a detector is the use of 
the well-known Fermi-Walker frame. If smeared detec¬ 
tors undergoing non-inertial trajectories are considered 
approximately rigid and sustained by internal cohesion 
forces, the smearing function becomes a function of the 
Fermi-Walker coordinates on the spacelike orthogonal hy¬ 
persurface to the detector’s motion. How this localization 
works out for general trajectories of the detector is very 
well detailed in |39j (or among other references, in, for 
example, [ISISilM])- 

Back to our original scenario, consider now two sta¬ 
tionary Unruh-DeWitt detectors A and B whose profiles 
are centred at positions Xf, and Xj^. They are mutually 
at rest and they share a proper time t. The interaction 
picture Hamiltonian describing the interaction of these 
two detectors with the field is 

Hi = '^XyXi^{t)fi^{t) f d^x F{x- x^)${x,t), (2) 

1/ 

where v = {A, B} labels the detectors, and the monopole 
moment operators act on the two-detector Hilbert 

space Hx <8) "Hb according to 

/lA(t) = TOA(t) 0 1 b, /iB(t) = 1 a 0 mB(t), (3) 

where 

m^{t) = ((T+e‘^"‘ -k (4) 

Let us now consider an arbitrary initial state density 
matrix po for the system comprised of the two detec¬ 
tors and the field. After the interaction of the detec¬ 
tors with the field, modulated by the switching func¬ 
tions Xy(t), the time-evolved state will be given by 
Pt = IJpoU\ where U is the interaction picture time 
evolution operator. Assuming that each is small 
enough to be in the perturbative regime , we can consider 
the Dyson expansion of the time-evolution operator U: 
U = U(o) -k [>(1) -k U(2) -k 0{Xl) where = 1 and 

poo poo pt 

dtH{t), dt dt'H{t)H{t') . 

J — oo ^ — oo ^ — oo 

( 5 ) 
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Hence pt = Po + Pt'* + Pt'^ + O(A^), where 

pW = c/(i)pp+PQC7(i)t ^ (6) 

pP = poU^^^^ + i7(2)po + PoU^^^^ . (7) 


The final state of the two-detector subsystem will, hence, 
be given by 


Pd,T = Tr0(pT) = Pd.o + Pd^T + Pd)r + : (8) 

where Tr^ denotes the trace over the held Hilbert space. 

As we will detail below, we are going to consider that 
detector A interacts with the held earlier in time, and 
detector B will couple to the held afterwards. We would 
like to analyze the ability of A to signal B. In particular, 
we would like to see what is the inhuence of the existence 
of A in the time evolved quantum state of B. Obviously, 
in a causal model, detector B’s state cannot depend in 
any way on the initial state of detector A if A and B 
remain spacelike separated during their interaction with 
the held. If we assume the most general uncorrelated ini¬ 
tial state for the detectors and the held, the initial state 
density matrix takes the general form po = Pd,o ® P^q^ 
where pd,o and p^ ^ are respectively the initial state of 
the two-detector subsystem and the initial state of the 
held. The leading-order contributions of the inhuence of 
detector A on detector B will be proportional to AaAb. 
The hrst order term in ([^ hence does not contribute to 
the signalling between A and B. Let us then focus on the 
second order contribution. Under the assumption that 
Pd — Pd,o ® P^Qi the second order contribution to the 
time evolved state of the two detectors will be given by 


/ OO nOO 

dt / dt' 

-OO J-oo 


/ OO pt 

dt / dt' Xip{t)Xv{t') 

-OO j —OO 

X fliy(t^fzYjit )^d,0 ^7 ^ ] 

/ OO pt 

dt / dt' Xip{t)Xvi^') 

-OO J —OO 

^ Pd, 0 p 7 ] (t (t) , t , 3?^, t] 


(9) 


Where W[x^, t, t'] denotes the pullback of the Wight- 
man function on the detectors’ smeared worldlines. 


terms in , that is, the terms that would vanish if one 
of the detectors is not present. With this in mind, we 
can decompose 0 as 

P'dr = '^AAePsigiiai -b ^ noise (H) 


^[2') 

where the terms pl, contain only contributions pro¬ 
portional to Ay, and as such they are local to each detec¬ 
tor and they are not involved in the flow of information 
from A to B. One quick way to see this is to notice that 
in the partial state of detector B, i.e. Pb,t = TrA/Sd.T, 
the only leading-order contributions that depend on the 
existence of detector A at all are proportional to AaAb. 

The terms pbg^^jj can be further simplihed assuming 
that the detector’s switching functions are compactly 
supported and that their supports do not overlap in time. 
Without loss of generality let us assume that the detec¬ 
tor A is switched on before B, and then let us make the 
additional assumption that 

supp[xA(i)] = [Tr^Tf], supp[xB(i)] = [r™,T°®], 

> Tf. (12) 

Imposing this assumption on we can simplify 
to 


/ OO poo 

dt / dt'XA{t)xB{t') (13) 

-OO J —OO 
X lU[a;B, t', a^A, t] 

+W[xt,, t, a;B, t'] {jj.B{t')pdfiPx{t)-pd,dPA{t)pLji{t') 


In a signalling scenario under the assumption (121, 
since A interacted with the held hrst, we can regard A 
as the sender and B as the receiver. All the informa¬ 
tion transmitted through the held from A to B will be 
encoded (to leading order) in the parts of the density ma¬ 
trix of B that are proportional to AaAb. This signalling 
contribution to B’s partial state will be given by 


J2) _ rp ,42) 

re,signal -^^Al^Pgignal/' 


(14) 


Recalling (|^ and making the reasonable assumption 
that prior to communication the initial state of the two 
detectors is uncorrelated (i.e., Pd,o = PaPb), we can 
compute the explicit form of (141 from (13): 


W[x^,t,Xn,t'] = J d^xj d'^x' F[x — Xiy] (10) 

X U[a:'-a:^]Tr^ [.^(a;, t)^{x', t')PAo] • 

Again Tr^ denotes the partial trace over the held Hilbert 
space. 

Since we are going to analyze the signalling between 
detectors A and B, we will not be interested in the local 


/ OO poo 

dt / dt'xA(t)XB(t')Tr[mA(t)pA] (15) 

-OO ^ —OO 


X 2ilm 


(^mB(t')PB-PBTOB(t') 


W[x^,t,xji,t'] 

where we have also used that 

lU[a:A, t, £Cb, t'] = {w[xji, t', Xj,, (16) 
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which is very easy to check form ( 101 . 


For each of the the individual detectors we may intro¬ 
duce a two-by-two matrix representation for each detec¬ 
tor’s individual Hilbert spaces. We will follow the same 
convention as in [49) . in which the ground |g) and excited 
|e) states of the two-level detectors correspond to 


and the monopole moments 


(17) 




0 


In this representation, the most general uncorrelated 
initial state for the two detectors is 


Pd,o — Pa Pb — 


O^A 


PI l-ftA 


“b Pb \ QQN 

p* 1 -a^y 


where a,y G R, P^ G C satisfying the conditions that make 
Pd,o a positive operator. One can now trivially evaluate 

Tr[mA(t)pA] = 2Re(/3Ae'^*‘) (20) 

and also the commutator 

[wb( ),Pb] _ 20,^) 2iIm(;3Be‘^“*') 

( 21 ) 


Substituting (20) and (21) in (15) we get 


/ OO POO 

dt / dt'xA(t)XB(t')R-e(/3Ae'^"‘) 

-OO ^ —OO 


X Im 




( 22 ) 


iOet'^ jg-iOet' 

-ie'^=*'(l - 2 aB) -21m(/3Be’^=*') 


2Im(/3Be’^=‘') (1 - 2aB) 


This expression can still be further simplified, since the 
imaginary part of the Wightman function can be ex¬ 
pressed in terms of the expectation value of the field 
commutator: 


lm(^Tr [(^(a;, t')P 0 .o]) 

^ (tt[${x, t)^{x', g] - Tt[p^ J{x', t')^{x, t)]) 
= 7Tr([(^(a;,t)(^(a;',t') - ^{x',t')p{x,t)]p^^^) 


1 

a 

1 

a 

1 


^{x,t),p{x',t') 


(23) 


Notice that this quantity is independent of the initial 
state of tl 
c-number. 


state of the field p^ q since the field commutator is a 


Now substituting (10) into (22) and using the identity 


(23) we can finally write 


= 2 

signal 


O POO 

dt / dt' XA{t)xB{P)Re {Pab"^^*) 

yo J —OO 


-21m(/3Be'*^=*') -ie-'^->*'(l - 2aB) 

igiOet'(1 _ 2a^) 2lm(^Be‘*^=*') 


(24) 


where C{t,t') is the final responsible for the causal be¬ 
haviour of the leading-order signalling, and it is simply 
the integral of the field commutator expectation over the 
detector’s smearing 


(18) C(t,t') = ij P^x'F[x — xP\F[x'— x^\ 


(25) 


${x,t),^{x',t') 


Notice that (24) is traceless, as it should be the case for 


all the different order perturbative contributions to the 
evolved density matrix of detector B. This stems from 
the fact that all the Dyson series perturbative correc¬ 
tions to the time evolved density matrix at given order 
in perturbation theory are traceless (see e.g. [ 12 ] for a 
proof). 

The full-density matrix leading-order expressions in 
Eq. (24) generalize for arbitrary smearings the results 
obtained in [5D| for the leading-order contribution to the 
transition probability of pointlike detectors. They also 
generalize the results obtained in m for pointlike de¬ 
tectors to smeared profiles. We notice already several 
interesting aspects of (24): First, all the leading-order 


contributions to the full density matrix of B coming from 
the presence of A are proportional to the smeared integral 
of the expectation of the field commutator, which is a c- 
number and, therefore, independent of the initial state of 
the field. The fact that the leading-order signalling con¬ 
tribution to the transition probability of a detector was 
independent of the background noise was already noted 
in |5D|. Here we show that this feature is not limited to 
transition probabilities and carries over to the full state 
of the detector, including the detector’s quantum coher¬ 
ences, and for arbitrary spatial smearing. Note, however, 
that this does not necessarily mean that it is equally easy 
to communicate (at leading order) using two-level detec¬ 
tors independently of the level of ‘noise’ in the field. This 
is not so because the noise terms in ( 11 ) (which give rise 
to contributions to Pb,t proportional to Ag) do indeed 
depend on the background field state. Hence, the signal- 
to-noise ratio will, in general, depend on the field back¬ 
ground state. 

Second, since the field commutator always vanishes for 
events that are spacelike separated, it is obvious from 


(24) and (25) that for compactly supported smearing 


functions, and if no UV cutoffs are introduced in the re¬ 
sponse of the detector, the time evolution of the full den¬ 
sity matrix of detector B is causal, as noted in a slightly 
different context in [36]. This is also an explicit leading- 
order extension for compactly supported smeared detec¬ 
tors of the results in all orders of perturbation theory in 
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m, which showed that, for pointlike detectors, the mi¬ 
crocausality of the scalar theory already guarantees that 
the flow of information from detector A to detector B 
is causal. In other words, here we see explicitly that if 
the two detectors are compactly smeared and spacelike 
separated during their entire interaction time, detector 
B is not sensitive to the state in which detector A was 
prepared or even to the very existence of detector A. 

Finally, note that (24) under the assumption (12) is 


better behaved in terms of divergences than the noise 
terms in (11). This is remarkable because it is well 


known that the second-order contributions to the time 
evolved density matrix present UV divergences for three 
spatial dimensions when the switching is discontinuous 
m- This divergence is present in the Ag noise terms, but 
it does not appear in the signalling contribution, even 
in 3-|-1-dimensional spacetimes and for pointlike detec¬ 
tors plus sudden switching [to see it one can check the 
expression of the field commutator in 3-1-1 spacetime di¬ 
mensions (A13) and how it enters in (24)]. We shall 


discuss this point further in section [III C[ where we will 
explicitly consider the pointlike limit of a scenario where 
both detectors are suddenly switched on and off in 3-1-1 
dimensions. 


III. EFFECT OF NON-COMPACT SMEARING 
ON THE CAUSALITY OF DETECTOR MODELS 

A. Explicit expressions of C{t,t') for Gaussian 
smearing 

We have discussed how the signalling between the two 
detectors is strictly zero for compactly supported space¬ 
like separated detectors due to the presence of the field 
commutator in (24). In this section we will see how 


one of the most common non-compact smearing profiles 
(namely Gaussian) affects the causality of the model in 
the interaction of two particle detectors with the field. 
Let us consider the following spatial smearing 


F{x) = 


1 


(cry^)" 


2 / 2 
/cr 


(26) 


where n is the number of spatial dimensions. Let us also 
recall the value of the field commutator in flat spacetime 
for 1, 2 and 3 spatial dimensions (calculated in full detail 
in appendix [A| : 

[(j){x, t), ^{x', t')] = ^sgn(t' - t)e (|t -t'\-\x-x'\) 

(27) 


for 1-1-1 dimensions, 

[${x,t),(f){x',t')] = 

i sgn(t' — t) 


for 2-1-1 dimensions, and 

[0(a;,t),(^(x',t')] = (29) 

I-^ W k(^ -t' + \x- x'\) - Sit -t' -\x- x'\) 

\x — x'\ An I 


for the 3-|-1-dimensional case (note that, as operators in 
the Hilbert space of the field, they are all multiples of the 
identity but we have not made this explicit to relax the 
already heavy notational load of this paper). 


With these ingredients we can evaluate (25) for the 


different flat spacetime dimensionalities. Let us consider 
that the two detectors are at rest and their centres of 
mass are located respectively at x^^ and a^g. Without 
loss of generality, we are going to assume that a^A = 0 
and therefore ja^gj = L is the distance between the two 
detectors’ centres of mass. 


1. Gaussian smearing in 1+1 dimensions 
For the 1-|-1-dimensional case, an explicit evaluation of 


(25) yields 


-| /'CO pco 

Ci(t,t')=-x - y/ dx (30) 

2TTa^J_^ 

X 0 (|t — t'\ — ja; — a;'|) 

we can make the following change of variables 

u = x + x', V = X — x' , (31) 

which makes the Heaviside function easier to handle. We 
evaluate this expression analytically recalling that, in our 
setting, detector B will be switched on after detector A 
[see (12) and (24)]. Taking this into account we obtain 


Ci(t,<') = 

_ 1 

“ 4 


_ pt' — t ^ _(u-v-2L)^ 

-^ du dve e (WP 

47rCTV-oo Jt-t' 

Erf( + ^ 


V2o 


V2o 


(32) 


2. Gaussian smearing in 2+1 dimensions 
For the 2-|-1-dimensional case 


€ 2 ( 1 , t') = - 


1 


2TT^a^ 


Id^x Id^x'e 


e[it- 1 ')^ - \x - x'\^] 


(33) 


27 r -^/it — t'Y — ja: — x'\'^ 


0[(t — t'Y ~ I® ~ we can perform a change of variables similar to (31) 


(28) 


U = X + X 


V = X — x' 


(34) 


































yielding 


which can be simplified to 


7 






au d^ve 77^7^ e 


X 


^(t-i.)2-|^|2 


(35) 


which can be simplified to 

r 2/_2 „ 

g-i 


C2{t,t') = - 


47r3(j4 


9 |'ur- 2 -ua;B / 9 | ■» |^+ 2 -u-a!B 

d ite d ve 




(36) 


^(i_i/)2_|„|2 

The integral over d^it can be readily solved analytically 

poo , , 0 /*27r , ,, 

/ 11 II I I'^r / l-^l^ pos y 

/ d|it||tA|e 2^ / cl(^e 
JO JO 

'L|m| 


= 2tt d|M||M|e 2 ^/g 


= 2707^62^ (37) 


where in the intermediate step loix) is the 0 -th modified 
Bessel function of the first kind. Substituting in (36) we 
get 


C2{t,t') = - 


e 2 o 


d|i;||u| 


e 2,7 


r 2 Tr 


dip e 


27r2CT27o - \v\^ Jo 

l 2 , \^\2 

e“ 2 ^ /■* “* d|t;||u|e“ 2 ^ f L\v\ 


7’'o-2 7o |u|2 V ^ 


which unfortunately, to the best of the author knowledge, 
does not admit a closed expression in terms of well-known 
functions. It is nevertheless regular and easy to evaluate 
numerically. 


3. Gaussian smearing in 3+1 dimensions 


Co{t,t') = - 


g r |„|2-2u.;bb f r, _ Mf+2ig^ 

d'^ue 2,72 d'^ve 2,72 


Stt^ctS 




(42) 


The integral over d'^it can readily be solved analytically 


27r J d|M||Mfe 2 ^ J d(cos 6 >)e ^ 

dTTCr^ 


m: 


d|M||M|e 2.72 sinh 


Lu 


= 2V^i 


cr 6 2^2 


(43) 


\v\L cos >p 


(38) 


substituting in (42) and performing the integration over 
d^u we get 


p~++ l•.7|2 

Coit, t') = d|u| \v\e-^S{t -t' + 1^1) 

-| 77 |I, COB e 

X J d(cos 6 >)e 

-e-i^ , _|7,!2 , JL\v\ 

d|u| e ^(5(f — t -I-|t)|) sinh 


y/JjJ^aLJo 


—e ^ _(t'-t)2 ^ — 

e 2,72 sinh ' 


V^TT^aL \ 

which can be further simplified to 


(44) 


-1 


[I.-(t'-t)]2 [L + (t'-t)p 


Co{t,t') = — - -— I e 2 

2V^<jL 


where the peak at the light cone is explicit. 


(45) 


For the 3-1-1-dimensional case we get that 


S{t — t' + \x — x'\) — S{t — t' — \x — x'\) 

(39) 


\x — X' 


we can again perform the change of variables (34): 
M = a; -I- a;b v — x — x' 


(40) 

we also have to disregard the contribution of one of the 
deltas recalling that, in our setting, detector B will be 
switched on after detector A (see (12) and (24)) so that 
in the whole integral evaluation of (24) onlyT' > t yields 
non-vanishing contributions: 


Co{t.t') = 




/• „ /• „ ++-V+ (77-7,-20,3)2 

Id'^u d'^ve (2‘')2 e <2^72 




B. Causality with a Delta switching 

It is possible to do a clean analytic study of the causal 
response of the Gaussian smeared detector model con¬ 
sidering detectors that interact with the field instanta¬ 
neously. Namely, we consider that detector A perturbs 
the field only at an instant and detector B switches its 
interaction at an instant This technically translate 
into assuming that the switching functions are given by 

XA(i) =<5(f-G), Xsit) = S{t - Ib) (46) 

where, according to the setup that we have been consider¬ 
ing tg > tA- With these switching functions, the leading- 
order signalling contribution from the presence of Alice 
to the time-evolved density matrix of Bob is given by 

/5Sig„ai = 2Re(/3.e'^*‘^)C(G,tB) 

-2Im(/?Be'^«*") -ie-'^=*“(l - 2aB) 

_ 2ag) 2Im(,dBe‘^”*°) 


X 


(41) 


X 


(47) 
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Figure 1. Estimator £^n(A, L, a) of the signalling contribution of the presence of detector A to the state of detector B for (from 
left to right) n = 1, 2 and 3 dimensions as a function of their spatial separation L (in units of the time separation between their 
interactions with the field, A) and for various characteristic leghtscales cr of the Gaussianly smeared detectors. The maximum 
light-contact between the two detectors occurs when the two detectors are fully null-separated L = A. The different lines 
correspond to values of cr = 0.02A (red triangles), a — 0.05A (green rhombi), a — O.IA (orange squares), cr = 0.2A (blue 
circles). We see that in the pointlike limit the estimator vanishes for spacelike separation between detectors. The timelike 
signalling is non-vanishing for timelike separation in 1-l-lD and 2+lD showing the violations of the strong Huygens principle. 


and, from (47), we can see that in this case the signalling 
contribution will in general be zero if and only if te) 
is zero. We discussed that if the detectors are compactly 
supported, this will happen when the detectors become 
spacelike separated, this is, if we dehne A = and 

recall that we defined L = \xji — Xj^\, for pointlike sup¬ 
ported detectors we have that 


A < L => te) = 0. 


(48) 


This is not the case for non-compact detectors for which 
never vanishes. To estimate how much casual- 
ity is violated in the communication between detectors 
A and B we build the following estimator of the causal 
influence of A on B in Eq. ( [47| ) 

= (49) 


to which all the matrix elements of detector B in (47) are 
proportional, modulo constants of order 1. From (321, 
(38), (45) we get that 


£:i(A,L,cr) = i 


Erf 


L-f A 


V2o 


-Erf 


L- A 

V2a 


for one spatial dimension, 

_ 9 

cA 


£ 2 (A, L,cr) = 


e 2 cr 


dy 


ye 


10 


a/A2 - J/2 

for two spatial dimensions and 


^3(A, L,cr) = 


2\/2TT^aL 




/o(4l) (51) 


{L+^y 


— e 


(50) 


(52) 


for three spatial dimensions. 

These magnitudes reflect the strength of the signalling 
from A in the full state of B after interaction, and thus 
they should be zero when A and B remain causally dis¬ 
connected, that is, when A < L. In Fig. [^we show, for 
different values of a, how the signalling estimator behaves 


as a function of the spatial separation of the detectors 
interacting for a hnite amount of time A. We see how 
in the pointlike limit cr —> 0, the model is fully causal 
and the signalling estimator vanishes when the detectors 
are spacelike separated, that is, when L > A. We ob¬ 
serve how for one and two dimensions there is a leakage 
of the ability of the detectors to communicate into the 
timelike area of the light-cone. This is due to the viola¬ 
tions of the strong Huygens principle [SUES]. This phe¬ 
nomenon enables the transmission of information without 
being supported by an energy flow between sender and 
receiver, as proved in [20]. Interestingly, although flat 
3-|-I-dimensional spacetimes do not allow for this kind of 
timelike energyless signalling, almost any other non-flat 
four-dimensional spacetime does, as for example in ex¬ 
panding universes |53|. For our purposes, we also see in 
Fig. [T] how for finite-sized detectors smeared in a non¬ 
compact way, the setup suffers from superluminal sig¬ 
nalling. We will quantify this violation of causality and 
see how fast considering a sufficiently small it is possible 
to recover approximate causality. 

This is done in the analysis performed in Fig. In 
the top panel of Fig. we see that for a given detec¬ 
tor spatial separation L = l.IA (where A is the spatial 
separation) the undesired faster-than-light signalling is 
reduced faster than exponentially as cr 0 regardless of 
the number of spatial dimensions. Thus, independently 
of the dimensionality of spacetime, it is possible to re¬ 
cover an approximately causal model for the quantum 
communication between two detectors when a L — A 
to arbitrary precision when considering non-compactly 
smeared detectors. 

On a further step, we analyze a crucially different case: 
in the bottom panel of Fig. we consider a scenario 
where two detectors are separated by a cr-dependent dis¬ 
tance. In particular we consider that L = A -|- 2cr. That 
is, the detectors are separated by a distance always equal 
to twice their characteristic size cr from the lightcone. If 
the detectors were compactly supported with size cr this 
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would mean the detectors would be spacelike separated 
but still close to each other in relation with their own 
size. The reason for this analysis is to understand if in 
order for the causality violations to stay in check it is 
also necessary that the spatial separation is bigger than 
the detectors’ characteristic length, that is, whether we 
need L ^ cr in order to recover approximate causality or 
not. 



cr/A 



(t/A 

Figure 2. (Top) Estimator £n(A, L, a) of the signalling con¬ 
tribution from detector A to the state of B when the cen¬ 
tres of mass of the two detectors are separated by a distance 
L = I.IA as a function of the smearing lengthscale a (in units 
of A) for n =1, 2 and 3 dimensions. We see that the signalling 
when (7 <C (b — A) (which would be zero for compactly sup¬ 
ported detectors) decays over-exponentially fast. (Bottom) 
£’„(A, L, a) for L — A + 2a: when the centres of mass of the 
two detectors are separated by twice the ‘size’ of the detectors 
from the spacetime point of maximum light-connection. For 
compactly supported detectors of size a it would imply space¬ 
like separation and thus £„ should be zero. We see 1) the over¬ 
exponential dampening of the signalling when cr —>■ 0, and 2) 
the signalling from A to B peaks at a dimension-dependent 
value of a before decaying polynomially on the lengthscale a 
(in units of A). The different lines correspond to the number 
of spatial dimensions: n = 1 (blue circles), n = 2 (orange 
squares), n = 3 (green rhombi). The inset shows a zoomed-in 
area for a/A < 0.8. 

We see in the bottom panel of Fig. that, indeed, 
causality is recovered faster than exponentially as a goes 
to zero also in this case, as expected from the analysis 
above. However, there is another aspect of the figure 


that results interesting and perhaps unexpected at first 
sight: the causality violations when L = A + 2a grow 
with the size of the detector only up to certain point, 
which depends on the number of spatial dimensions, be¬ 
fore starting a slow polynomial decay as a keeps growing. 
This perhaps surprising effect stems from the fact that 
there is a distance decay in the integral of the commuta¬ 
tor in the spacelike area of the field in all cases (explicit 
in 2-1-1 and 3-1-1 dimensions but also existent in the 1-1-1- 
dimensional case) This decay is not dependent at all on 
the detector size lengthscale a but it does depend on 
the dimensionality of spacetime. This means that, for 
extremely large non-compact supported detectors whose 
centres of mass are spacelike separated by only twice their 
characteristic distance, it is possible to recover some de¬ 
gree of approximate causal behaviour in principle, albeit 
only very slowly converging to causality as the size of the 
detectors increases. 


C. Dependence on the switching: The regular 
behaviour of the signalling contributions 


In the section above we have considered a delta switch¬ 
ing corresponding to the instantaneous interaction of the 
detectors with the field to keep the analysis clean. How¬ 
ever, one may legitimately wonder how considering dif¬ 
ferent compactly supported switching functions (includ¬ 
ing some finite-time interaction of the detectors with the 
field) may modify these results about the causality of the 
Gaussianly smeared detectors. 

To convince ourselves that the shape of the switching 
is not relevant for the results presented in this section, we 
can analyze the case of an extended switching function 
in 3-1-1 dimensions which again yields closed expressions 
for our signalling estimators. Furthermore, this analysis 
will allow us to show some interesting aspects about the 
UV regularity of the signalling contribution (p4|) . 


For simplicity, let us assume for this subsection that 
the detectors are gapless. This assumption will remove 
from the study undesired oscillations coming from the in¬ 
ternal detector dynamics that are irrelevant for the causal 


behaviour of the detectors. From (24) we can see that 


a reasonable signalling estimator for compact switching 
functions with non-overlapping support can be 


fc = 


dt / dt' XA(t)xB(t')C(t,t') 


(53) 


Let us consider the following simple sudden switching 
scenario where the interaction of both detectors has a 
duration T: 


xAt) 


1 t G [t^, ty T\ 

0 t (fi tu -\- T] 


(54) 


Without loss of generality we consider the switch-on time 
of detector A = 0 and we call A the time interval 
between the switch-off of detector A and the switch-on 





















10 


of detector B. In that case we can express (55) for the 


3+1-dimensional case, using ( |45| ), as 
1 


fc = 


dt / dt e — e 2 ^^ 


2y/2'K^aLJo Ja+t V 
which admits the following closed form: 

£c(A,L..,r, = j+[(L-A)Erf(^ 

- (A+L)Erf ( ) +2(A-L+r)Erf 


V2o 


(55) 

(56) 

L-A-T 

V2a 


+ {L- A- 2T)Erf (- - "" ) + 2(A + L + T) 


'L- A-2T 


\ 


'A- 




vIk 




(A±^ 


(A-.L + T)^ 


— cr e —e 2 a^ —e 


(A + -L + T)-^ 

9^2 


(A-L+T)^ 


— e 


+e 


(A-.L + 2T)^ _ (A + -L + T) 

2 ^^ -l-e 2 ^^ 


^ +e 

^ _ (A + -L + 2T)^ 

— e 2 ^^ 


In Fig. I^we see how the behaviour is qualitatively the 
same as in the case of a delta switching just accounting 
for the larger amount of causal contact due to the du¬ 
ration of the interaction T. In the figure, the signalling 
contribution is smaller in magnitude than in the instan¬ 
taneous switching scenario (which causes a more violent 
perturbation of the field). Notice an important point: as 
it can be seen in the plot, in Fig. the pointlike limit 
CT —?► 0 is no longer a delta distribution, but instead a 
finite function that has support from the point when the 
causal contact is initiated at L = A, peaks at maximum 
causal contact L = A + T and decreases down to zero 
when the causal contact finishes at L = A + 2T. 

We can actually compute the exact limit cr —> 0 on 


(56), which yields 


lim £c{A,L,(7,T) 

cr—>0 


(57) 


|L-A| + |L-2T-A| -2|L-T-A| 
47rL 


The signalling estimator peaks at maximum light-contact 
L = A + T yielding a finite value, even in the pointlike 
case, of 




max 


1 

27r(I +A/T)’ 


(58) 


which also does not diverge when the duration of the 
interaction is considered large T —> 00 . 

This is a very remarkable aspect of the leading-order 
signalling contribution from a detector A to a detector 
B: It does not diverge when a non-continuous switching 
function and pointlike detectors are considered. This is 
in stark contrast to the known divergences in the noise 
terms proportional to Ag for 3+I-dimensional spacetimes 
m- What is more, this signalling contribution remains 


finite and controlled even for long interaction times. This 
means that, in what concerns the leading-order terms, 
the signalling contribution can be perturbatively investi¬ 
gated in a consistent way even for very long interaction 
times. This is however only fully consistent if the larger 
order corrections to higher power signalling terms were 
also to remain finite in the limit T —)■ 00 . Regarding the 
causal behaviour of the model, we see that the acausal 
predictions of the non-compact smeared detectors with 
a finite sharp switching behave in an identical way as 
the case of the delta switching, reflecting the switching 
independence of the observed behaviour. 



L/A 

Figure 3. Estimator £c{A, L, a, T) of the signalling contribu¬ 
tion of the presence of detector A to the state of detector B in 
a 3+lD spacetime when the detectors are sharply switched on 
and off after allowing them to interact with the field for a finite 
time (T = 0.1 A) as a function of their spatial separation L (in 
units of A) and for various characteristic lenghtscales u of the 
Gaussianly smeared detectors. The maximum light-contact 
between the two detectors occurs when L — A + T . The 
different lines correspond to values of cr = 0.02A (red trian¬ 
gles), cr = 0.05A (green rhombi), a = O.IA (orange squares), 
a = 0.2A (blue circles). In the pointlike limit the estimator 
vanishes for spacelike separation between detectors. The sig¬ 
nalling contribution is finite in this limit even with a sharp 
switching. 


IV. EFFECT OF A UV-CUTOFF ON THE 
CAUSALITY OF COMMUNICATION 

In this section we will analyze the violations of causal¬ 
ity in the communication between two particle detectors 
steaming from the introduction of a hard UV cutoff A 
in the interaction between the field and the detector in 
the privileged quantization frame (t,x). As discussed 
in the introduction, it is common to carry this kind of 
approximations when modelling the light-matter inter¬ 
action. The rationale of such an approximation is the 
fact that the coupling between an atom and the elec¬ 
tromagnetic field is frequency dependent and will even¬ 
tually vanish for large enough electromagnetic energies 
(the atoms become transparent to all frequencies above 
certain threshold). 
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This approximation, of course, strongly affects the 
causal behaviour of the model. In this section we will 
perform a clean study of the emergent causality viola¬ 
tions that appear in the two-detectors communication 
model when a hard cutoff is introduced. Since we have 
to have a clean signature of the effect of the UV cutoff 
and not mix it with the effect of a spatial smearing, we 
are going to consider delta switching and pointlike de¬ 
tectors for this study in the cases of H-l, 2-1-1 and 3-1-1 
dimensions. 


The hard UV cutoff is implemented by directly sev¬ 
ering the interaction of the field with the detector for 
field frequencies above a given threshold frequency A. 
To model this, we need to replace the field operator 
that couples to the atoms in Q by an operator 
which incorporates an UV cutoff A in its mode expansion. 
Regarding the detector’s causal response, this cutoff in¬ 
volves substituting the held commutators in (24) with 
UV cutoff versions. 


1+1 dimensions 


For the case of one spatial dimension we start from 


(A4) and apply the UV cutoff: 


= ^ ( / d/c sin [k{t — t' — x + x')] 

nO ^ 

+ dk —- sin [k{—t + t' 
J-A 


(59) 



The integrals in the expression above are exactly the def¬ 
inition of the sine integral functions, so we can rewrite 
this UV-cutoff commutator as 


[$a{x, t),^A{x',t')] (^Si[A(t -t' + x- x')] 

— Si[A(< — — X + a;^)]^ (60) 


where Si(z) is the sine integral function. 

As mentioned above, we are going to assume that the 
detectors are pointlike and placed at positions Xa, Xb, and 
that they interact instantaneously with the held at times 
U and ts- Namely, 


XAit) = S{t - tj,), Xsit) = S(t - fe), F{x) = S{x) (61) 


in Eq. ( |24[ ). Without loss of generality, let us assume that 
Xa — 0,Xb = L, Ia = 0 and = A. With these choices, 
we can write the signalling estimator (49) dehned in the 
previous section as 


1 


£-f(A,L) = —|Si[A(L + A)] -Si[A(L-A)]|. (62) 


2+1 dimensions 


We operate under the same hypotheses as in the 1+1- 
dimensional case (pointlike detectors and delta switch¬ 
ings). Upon the introduction of a UV cutoh in the ex¬ 
pression for the held commutator (A8) and following the 


same prescriptions as above we can compute the form 


of the signalling estimator (49) for the 2+1-dimensional 
case 




rA 


dy Jo {yL) sin [y A] 


(63) 


which cannot be evaluated i a closed form, but which is 
again easy to evaluate numerically. 


3+1 dimensions 


We repeat the same analysis for the 3+1-dimensional 
case. From (A12), the cutoh version of the held commu¬ 
tator is 


[ij)A[x,t),ij)A{x' ,t')\ = 


Stt^Ix — x'\ \ J_A 
A 




-J (64) 


which can be readily evaluated to the following closed 
expression 


^sin [A{t — t' + \x — x'D] sin [A(t'— t + \x — x' 


(65) 


t — t' + \x — x' 


t' — t + \x — x' 


Again let us call L = \xa — Xb\, and A = te — ^a- With 


these choices, we can write the signalling estimator (49) 
for the 3+1-dimensional case as 


S^{A,L) = 


1 




sin [A(L — A)] sin [A(L + A)] 


L-A 


L + A 


( 66 ) 


We show in Fig. the behaviour of the estimator 
£^{A,L) as a function of the spatial separation of the 
detectors L in units of their time separation A for dif¬ 
ferent increasing values of the UV cutoh A and for 1+1, 
2+1 and 3+1 dimensions. We see in which way the fully 
causal model is approximated as the cutoh increases. 

One can quickly appreciate a rather concerning as¬ 
pect of the causal behaviour of this model in the three- 
dimensional (space) case: we can see from the analytic 
expression (66) that increasing the cutoh A will not fix 


the acausal signaling for a fixed distance L when consid¬ 
ering spacelike separated detectors {L > A). Instead of 
restoring causality, increasing the cutoh will induce faster 
oscillations, but the magnitude of the acausal influence 
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Figure 4. Estimator S^{A, L) of the signalling contribution of the presence of detector A to the state of detector B for (from 
left to right) n =1,2 and 3 dimensions as a function of their spatial separation L (in units of the time separation between their 
interactions with the field, A) and for various UV cutoffs A. The maximum light-contact between the two detectors occurs 
when the two detectors are null-separated L = A. The different lines of the 1-l-lD plot correspond to values of A = SA”'^ 
(blue circles), A = 10A“^ (orange squares), A = 20A“^ (green rhombi) and A = 10®A“^ (orange triangles). For the 2+lD 
and 3-l-lD plots they correspond to values of A = 25A“^ (blue circles), A = 50A“^ (orange squares) and A = 100A“^ (green 
rhombi). Since the detectors are pointlike and interact with the field instantaneously, the non-vanishing values in the spacelike 
separation area {L > A) are due to causality violations caused by the hnite UV-cutoff. Signalling is non-vanishing for timelike 
separation in 1-l-lD and 2-l-lD even in the limit A —^ oo due the violations of the strong Huygens principle. 



AA 


Figure 5. Scatter plot of the estimator £^[A,L) of the sig¬ 
nalling contribution of the presence of detector A to the state 
of detector B for a fixed spacelike separation of the detectors 
L = 1.4A as a function of the UV cutoff scale A (in units 
of A~^) for spatial dimensions equal to n = 1 (blue circles), 
n = 2 (orange squares), n = 3 (green rhombi). In one and 
two spatial dimensions causality is restored polynomially in 
A. In three-dimensional spaces increasing the cutoff does not 
restore causality unless we go to the strict limit A —>■ oo. 

of A on detector B will remain unaltered and causality 
will only be recovered in the strict limit A —>■ oo. This 
can be seen graphically in Fig. where we see that in¬ 
creasing the cutoff reduces the acausal influence of A on 
B at a fixed L polynomially with the UV cutoff for 2-1-1 
and l-bl dimensions, but it does not reduce it for the 
3-1-1-dimensional case. 

In particular, for 2-1-1 and 1-1-1 dimensions we see that 
the acausal influence decreases as ~ A““, with a > 1 for 
1-1-1 dimensions and a < 1 for 2-1-1 dimensions. These 
results are consistent with the result for 1-|-1-dimensional 
Dirichlet cavities obtained in [19] where the causal influ¬ 
ence of one detector on another decayed with the square 
of the number of cavity modes considered in the model. 


Although the ill behaviour in 3-1-1 dimensions may 
of course be, to some extent, smoothed by considering 
smoother switchings than a delta-kick, the fact remains 
that the the restoration of the casual behaviour of the 
model in 3-1-1-dimensional spacetimes is much worse than 
in 1-1-1 and 2-1-2 dimensions, and so the use of hard cut¬ 
off detector models in three spatial dimensions for detec¬ 
tors switched on for a finite time to scenarios where the 
causality of the model is relevant is problematic. 

We conclude that while it may be argued that approxi¬ 
mate causality may be recovered with arbitrary precision 
polynomially fast in 1-1-1- and 2-|-1-dimensional scenarios 
if we introduce a hard UV cutoff, it is not a reasonable as¬ 
sumption to introduce UV cutoffs in relativistic quantum 
communication scenarios in 3-1-1-dimensional spacetimes. 

V. CAUSALITY VIOLATIONS OF THE 
ROTATING-WAVE APPROXIMATION 

The rotating-wave approximation (RWA) is yet an¬ 
other very common approximation made in the modelling 
of quantum optical settings. In fact, it is arguably the 
most common approximation in quantum optics, and it 
can be ubiquitously found anywhere from basic textbooks 
to research works [T^j. To better understand this approx¬ 
imation, let us first expand the held in plane-wave modes 
in the Hamiltonian Q: 

Hi = j (67) 

where the integral over d"a; has already been performed 
and where we have dehned the Fourier transform of the 
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real-valued spatial profile as 
F{k) = ^ 




d"a: F{x)e 


ik-\ 


( 68 ) 


The RWA consists of removing the terms proportional 
to and their Hermitian conjugates in the Hamil¬ 

tonian (67). More explicitly, the so-called rotating-wave 


approximation consists of replacing the Hamiltonian (67) 

by 




Hk) 

V\k\ 


(69) 

The rationale behind this approximation is that the ne¬ 
glected terms do not have a stationary phase for the 
detector-held resonance fc = H, and as such, those 
terms yield bounded oscillations when integrated in time. 
Within the approximation, these bounded oscillations 
can be neglected in the detector-held dynamics when 
compared to the close-to-resonance rotating-wave terms. 
Thus, the Hamiltonian (69) is expected to approximate 


the light-matter interaction for interaction times longer 
than the Heisenberg time of the detectors -1 m- 

Notice that from the form of this Hamiltonian we al¬ 
ready see that if the initial state of the held is the vacuum 
and the detectors are both in the ground state, there is 
no non-trivial dynamics: the detectors will not be able 
to communicate or get excited even if the interaction is 
switched on sharply and for a hnite time. 

Since the Hamiltonian is no longer linear in the held, 
the microcausality of the scalar held theory no longer 
guarantees that the model will behave causally at all. 
To illustrate and study the causality violations in this 
approximated model, we are going to consider a partic¬ 
ular scenario where we begin with a general state of the 
two detectors, A and B, but the held is initialized in the 
vacuum. This is akin to the classical scenario of com¬ 
munication with a pair of antennas in the presence of no 
background noise, 


/3o=/5d.o®|0)(0|. (70) 


To have a clean result and disregard all the other pos¬ 
sible sources of acausal behaviour that we have already 
studied, we are also going to consider pointlike detectors. 

The leading-order signalling will have, again, two con¬ 
tributions corresponding to each term in ([^ that is pro¬ 
portional to AaAb. We will focus on the 3-1-1-dimensional 

case and let us evaluate Tr^ ): 


Tr^ E ( 71 ) 


IX, 

3/ 


xy df y dt'xAt)Xvit')J 


where we have used that Tr^ (al |0)(0| Uk') = S^^^(k — k') 
and the fact that the detectors are considered pointlike, 
that is F(k) = 1. 

We can quickly evaluate the integral over d^fc that we 
will denote as 


rft f' T x ) — k [\(^\k\-n){t-t')-k-(x^-Xr,)] 

j 1^1 c. 




•^771 JO 


d|fe| ) sin (|fc||a;i, — a:,j|) 




ITT 

-b — 


— Xril"^ — {t — t'Y 2 
-t' - \Xy - Xr^\) - 6{t -t' - \x^ -b Xj^\) 


(72) 


Here we see explicitly how for the signalling terms (pro¬ 
portional to AaAb), the rotating-wave approximation 
breaks causality: if we compare this expression with 


the commutator (29), we see that the result under the 


rotating-wave approximation is essentially the same in 
terms of causal behaviour than the full model without ap¬ 
proximation, except for the hrst summand in (72), which 


is explicitly non-vanishing outside the light-cone and thus 
introduces causality violations. 

Let us compute the signalling contribution from A to 
the density matrix of B in the terms Tr^ 
of ^ under the RWA by tracing out detector A in the 
terms proportional to AaAb of (72): 


Tr0,A = ^^^Tr(CT^/5A.o)/5B,oO-^ 

/ CXD pOO 

dt / dt'xj,{t)xB{t')C{t,t\xj^,Xji) 

-00 J —00 

+ Tr(pA,od+)(j“/5B,o (73) 

/ OO pOO \ 

dt / dt'xB{t)xAit')C{t,t',Xs,Xj,)Y 
-00 J —00 / 

It is easy to see that the signalling contribution to 
the density matrix of B coming from + Pol7(2)t 

yields identically behaved terms in terms of the scaling 
of causality violation as those coming from p^U^^d 
under the hypothesis (12). 


Therefore, to estimate the causality violations in the 
two-detector communication scenario coming from the 
rotating-wave approximation, we need to assess the mag¬ 
nitude of the following signalling estimator which will be 
zero if there is no influence of A on B: 


^RW — 


dt / dt'x,,{t)xB{t') 




L'^-it- t'Y 


(74) 


This is, the terms of the causal influence of A on B (Eq. 
(HH)) that are non-zero outside of the lightcone. Because 
it is argued that the rotating-wave approximation only 
becomes accurate for long durations of the interaction. 
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Figure 6. (Top) Estimator £rw of the signalling contribu¬ 
tion from detector A to the state of B under the RWA for 
two pointlike detectors separated by a distance L (in units 
of The duration of the interaction for both detectors 

is T = 17“^ and the interval between A’s switching off and 
B’s swithcing on is A = 17“^. The values of L where the de¬ 
tectors’ interactions have some light contact are enclosed by 
two vertical red dashed lines. The signalling estimator peaks 
at maximum causal contact. We see that the signalling at 
strictly spacelike separation {L > A + 2T) is non-zero due 
to the causality violation introduced by the RWA, decreasing 
slowly with the spatial separation. (Bottom) Srw for fixed 
L = 3T + A (detectors fully spacelike separated) as a func¬ 
tion of 12 in units of T~^: We see that while it is true that 
causality gets restored when T ^ 17“^, the causality violation 
is only polynomially recovered when 17 grows. 


which maps the integral to 

p2T+A pu-2T-2A 

Srw =7: du dve-^^^ 

I Jt+a J-u 

^ p — u-\-‘2T ^ 

+ du 

J 2 T+A Ju-iT-2A 

The closed expression is lengthy (a long sum of sine inte¬ 
gral and cosine integral functions) but easily evaluable. 

The two pointlike detectors are completely spacelike 
separated for values of L satisfying L > 2T + A. To 
see how the violation of causality of the rotating-wave 
approximation behaves with the increase in the duration 
of the interaction, we are going to fix L = 3T -|- A so that 
the detectors are always fully spacelike separated during 
the interaction. 

We see in Fig. |^that when T ^ 17“^ (the interaction 
times where the RWA is usually considered to start be¬ 
ing unacceptable), the causality violation decreases very 
slowly with the spatial separation betweeen the detec¬ 
tors. The fact that the RWA model allows for long- 
range faster-than-light signalling renders the model in 
this regime completely unsafe for any kind of relativistic 
considerations. 

We see also that as we increase the product flT the 
RWA violation of causality becomes smaller as expected. 
However, importantly, as flT increases, the speed at 
which causality is restored in the model becomes increas¬ 
ingly slower: while it is possible to improve the causality 
of the RWA-approximated model if the interaction time 
is much larger than the speed at which causality is 
recovered drops fast as the duration of the interaction in¬ 
creases. While it seems possible to work under the RWA 
and have some approximate degree of causal behaviour 
for very long interaction times, the spacelike signalling 
that the RWA enables is resilient, long-ranged (in terms 
of the spatial separation of the detectors) and remains 
even for relatively long timescales. This means that one 
should be extremely careful when using this approxima¬ 
tion in settings where faster-than-light signalling cannot 
be tolerated. For instance to make predictions about ex¬ 
periments where the outcome relies on the fact that the 
detectors remain spacelike separated. 




we are going to consider a finite-duration interaction. 
Namely, we will use the switching function ( [5^ . In other 
words, the interactions of A and B have a duration T, A 
switches on at t = 0 and there is a delay A between the 
switching off of A and the switching on of B. This leaves 
([74|) as follows 


Srw = 


dt 


I.2T+A inp-t') 


dS 


T+A L'^ — {t — t'Y 


(75) 


this integral can be evaluated analytically, for which we 
first perform the change of variables 

u = t + t', V = t — t\ (76) 


VI. CONCLUSIONS 

We have analyzed, in an operational way, the causal 
behaviour of three of the common approximations and 
generalizations of the particle detector models used in 
quantum field theory to probe quantum fields (such as 
the Unruh-DeWitt detector) and of the closely related 
models of light-matter interaction (such as the Jaynes- 
Cummings model). Namely, we have studied 1) the in¬ 
troduction of a non-compact spatial smearing of the de¬ 
tector, 2) the introduction of a UV cutoff in the detector- 
field interaction and 3) the use of the rotating-wave ap¬ 
proximation (RWA). 
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Concretely, we have studied how the three modifica¬ 
tions of the model enable unacceptable faster-than-light 
signalling between two particle detectors interacting with 
a quantum field at different times in fiat spacetimes of 
one, two and three spatial dimensions. While all these 
three modifications of the standard UDW model intro¬ 
duce some degree of pathological faster-than-light sig¬ 
nalling in quantum communication scenarios, we have 
seen that some are more harmful than others. 

First, we have analyzed the causality violations that 
appear when particle detectors are spatially smeared with 
a smearing function that is not compactly supported. 
This is of interest because this kind of approximation is 
very often considered in the literature, both as a means 
to regularize divergences or as a refinement of the light- 
matter interaction models to include the finite size of 
the (non-compactly supported) atomic wavefunction. We 
have seen that the causality violations introduced by a 
Gaussianly supported spatial profile can be made arbi¬ 
trarily small over-exponentially fast as we reduce the size 
of the detectors in units of the time separation between 
them, regardless of the number of spacetime dimensions. 

Second, we have studied the effect of the introduction 
of a UV cutoff on the signalling between two spacelike 
separated detectors. This kind of approximation is made 
in the context of the light-matter interaction and it is 
often justified by the fact that atomic probes do not cou¬ 
ple with the same strength to all the frequencies of the 
electromagnetic field, being effectively transparent to fre¬ 
quencies much above the atomic frequency resonances. 
We have seen that the effect of this kind of cutoff is crit¬ 
ically dependent on the dimension of spacetime, and so 
are the limitations on the model that relativistic causal¬ 
ity imposes. Namely, while the violation of causality can 
be made arbitrarily small for 1-1-1 and 2-1-1 dimensions 
polynomially fast on the frequency of the cutoff, the case 
of 3-1-1 dimensions is different: for very fast switchings it 
is not possible to make the faster-than-light signalling ar¬ 
bitrarily small for any finite value of the UV cutoff. This 
renders this kind of approximation in 3-1-1 dimensions 
unsuitable to describe scenarios in relativistic quantum 
communication or any other scenarios where relativistic 
considerations are important. 

Finally, we have assessed the violations of causality 
that emerge from the application of the rotating-wave 
approximation (commonly used in the light-matter in¬ 
teraction models such as the Jaynes-Cummings model). 
We have seen that the RWA allows for long-range faster- 
than-light signalling for interaction times of the order of 
the detector’s Heisenberg time The violations of 

causality do get reduced as the interaction time becomes 
larger than the detector’s Heisenberg time, but their re¬ 
duction is slow and that makes the RWA perilous to use 
and arguably an approximation to be avoided —even for 
long interaction times— when describing setups in rela¬ 
tivistic quantum communication and any other circum¬ 
stance where the causality of the model is crucial. 

Additionally, we have discussed some remarkable as¬ 


pects of the quantum communication between two parti¬ 
cle detectors through a quantum field. First, the leading- 
order signalling contributions of one particle detector to 
the state of another detector (that is, the influence that a 
detector A switched on at early times has over a detector 
B switched in A’s future) is completely independent of 
the background initial state of the field regardless of the 
initial state of detectors A and B. This is commensurate 
with what was seen for the transition rate of detector B 
in [20]. Nevertheless, we have discussed that signalling 
from A to B is still affected by the background state 
through the local noise terms. Second, we have seen that 
the influence of detector A on the state of detector B (re¬ 
ferred to as the signalling contributions throughout this 
work) does not fundamentally present the kind of UV 
divergences described in |50j for pointlike detectors and 
discontinuous switching functions. 
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Appendix A: Field commutators in different 
spacetimes 

In this appendix we obtain in full detail the scalar 
field commutator in fiat spacetime for H-l, 2-1-1 and 3-1-1 
spacetime dimensions. Although very basic, these calcu¬ 
lations are going to be useful to show how to obtain the 
field commutator when an UV cutoff is introduced. 

Notice that throughout this appendix the results of 
these commutator evaluations have to be regarded as c- 
numbers (multiples of the identity in the field’s Hilbert 
space). 


1. Free scalar field commutator in 1 + 1-dimensional 
Minkowski spacetime 


The easiest way to obtain this commutator is to con¬ 
sider a plane-wave mode expansion of the field: 

^{x,t)= [ dk / . 

J-oo \/47r|fc| ^ 


(Al) 


Using this we get 
[(j){x,t),}{x',t')] = 

f dk f dh' ^ fin, ]p-di\k\t-\k'\t')-kx-k'x')] 


1 

47r 




(A2) 
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where we have used that ak,ak' = =0. Now 


ak,al,] =5{k-k'), [a|,,afc'] =-[a\,,ak\ =-S{k'-k) 


we have that 


where we have used that [afc,afc'] = [ 0 ^, 0 ^/] = 0. Now 
since [afc,aj^,] = - k'), 


[^{x, t), Ux', t')] = ^ J^dk ^ I ' 


_ gi[(|fc|(*-d)-fe'(!®-^')] 


— i[(|fc| (i ——|fc| |£C —£C^ I COS (^] 


i[\k\{t-t')-k(x-x’)] 


= —f /" dfc ^ sin [k{t — t' — X + x') 
2tt\Jq k 


+ J dk sin [k{—t + t' — x + a;')] ^ 


f°° I TT 

Now since / dx—sin(ax) = —sgn(a), 

Jo X 2 

[^{x,t),lj){x\t')\ = y (^sgn[t-t'-x + x'] 

— sgn \t' — t — X + x'^ (A5) 

which can be simplified to 


Jj){x, t), J){x', t')] = ^sgn(t' - t)e (|t - t'l - |x - x'l) 


2. Free scalar field commutator in 2+1-dimensional 
Minkowski spacetime 


Again we make a plane-wave mode expansion of the 
field, as in the 1+1-dimensional case: 


^{x,t) = 


'{2TTy2\k\ 


We can now use this expansion to obtain that 


${x,t),(j){x\t') = 


_ \ cost/?] 


1 

= ^y d\k\Jo{\k\\x-x'\)sin[\k\{t-t')] (A9) 


which can be evaluated in closed form as 


[${x,t),${x',t')] = 

i sgn(t' — t) 

27r y^{t — t'Y — I® — 


= Q[{t-t'f -\x-x'\^ 


3. Free scalar field commutator in 3+1-dimensional 
Minkowski spacetime 

Once again, we expand the field in terms of plane-wave 
modes as we did in the previous cases: 


^(27r)32|fe| 


(All) 

With this expansion at hand we obtain that 
[J){x,t),${x',t')] = 

-^[d^k fdJk' , ^ 

leW J 

where we have used that [afc,afc'] = [afc,aj[,,] = 0. Now 
since [afc,aj^,] = S^^^k - k'), 


(A7) 








dJk / d^fc'- 


dfc, dk' ed(\H^-W\t')-k-x-k'-x')] 


+ [Ofc, dl,] e-i[(|fe|‘-|fe'|i')-fe-a!-fc'-a3')] 


= ^ y sin [\k\(t -t')-k-{x- x')] 

= Jd\k\\k\J d(cosd) sin [\k\[{t — t') — \x — x'\ cosd]] 

= i/d|fc|^^(cos[|fc|(t-t' + |x-a^'|)] 

— cos [\k\{t — t' — \x — x'D] ^ 
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i 1 


\x — x'\ Stt^ 


d|fc| 


+ \x-x' \) 


+ 


d|fc| 




d|fe| 


i|fc|(t-t'-|x-x'|) 


i 1 
la: — a:'| 


Jo 

dk 

D 


which finally yields the following closed expression with 
support only on the null boundary of the lightcone, (in 
contrast to the timelike leakage of 1+lD and 2+lD) 


[${x,t),^{x',t')] = 
i 1 


la; — x'l 47r . 


(A13) 

S{t — t' + |a; — a;'|) — S{t — t' — \x — x' 


[1] A. Dragan, J. Doukas, and E. Martm-Martmez, Phys. 
Rev. A 87, 052326 (2013) 

[2] A. Dragan, J. Doukas, E. Martm-Martmez, and D. E. 
Bruschi, Classical and Quantum Gravity 30, 235006 
(2013) 

[3] S.-Y. Lin, Ann. Phys. 351, 773 (2014) 

[4] W. G. Unruh, Phys. Rev. D 14, 870 (1976) 

[5] B. DeWitt, in General Relativity: An Einstein Centenary 
Survey, edited by S. W. Hawking and W. Israel (Cam¬ 
bridge University Press, Cambridge, 1979). 

[6] S.-Y. Lin and B. L. Hu, Phys. Rev. D 76, 064008 (2007) 

[7] E. G. Brown, E. Martm-Martmez, N. G. Menicucci, and 
R. B. Mann, Phys. Rev. D 87, 084062 (2013) 

[8] D. E. Bruschi, A. R. Lee, and 1. Fuentes, J. Phys. A: 
Math. Theor. 46, 165303 (2013). 

[9] L. G. B. Grispino, A. Higuchi, and G. E. A. Matsas, Rev. 
Mod. Phys. 80, 787 (2008). 

[10] S. Takagi, Prog. Theor. Phys. Suppl. 88, 1 (1986). 

[11] P. Candelas and D. W. Sciama, Phys. Rev. Lett. 38, 1372 
(1977) 

[12] M. O. Scully and M. S. Zubairy, Quantum Optics (Cam¬ 
bridge University Press, 1997). 

[13] A. Wallraff, D. 1. Schuster, A. Blais, L. Frunzio, R. S. 
Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. 
Schoelkopf, Nature 431, 162 (2004) 

[14] E. Martm-Martmez, M. Montero, and M. del Rey, Phys. 
Rev. D 87, 064038 (2013) 

[15] A. M. Alhambra, A. Kempf, and E. Martm-Martmez, 
Phys. Rev. A 89, 033835 (2014) 

[16] E. Martm-Martmez, D. Aasen, and A. Kempf, Phys. 
Rev. Lett. 110, 160501 (2013) 

[17] E. Martm-Martmez and C. Sutherland, Phys. Lett. B 
739,74 (2014) 

[18] M. Gliche and A. Kempf, Phys. Rev. A 81, 012330 (2010) 

[19] R. H. Jonsson, E. Martm-Martmez, and A. Kempf, Phys. 
Rev. A 89, 022330 (2014) 

[20] R. H. Jonsson, E. Martm-Martmez, and A. Kempf, Phys. 
Rev. Lett. 114, 110505 (2015) 

[21] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 
2738 (1977) 

[22] L. J. Garay, M. Martm-Benito, and E. Martm-Martmez, 
Phys. Rev. D 89, 043510 (2014) 

[23] E. Martm-Martmez and N. G. Menicucci, Glassical and 
Quantum Gravity 29, 224003 (2012) 

[24] A. Valentini, Physics Letters A 153, 321 (1991) 

[25] B. Reznik, Foundations of Physics 33, 167 (2003) 


[26] B. Reznik, A. Retzker, and J. Silman, Phys. Rev. A 71, 
042104 (2005) 

[27] G. Ver Steeg and N. G. Menicucci, Phys. Rev. D 79, 
044027 (2009) 

]28] S. J. Olson and T. C. Ralph, Phys. Rev. Lett. 106, 
110404 (2011) 

]29] G. Salton, R. B. Mann, and N. G. Menicucci, New J. 
Phys. 17, 035001 (2015) 

]30] A. Pozas-Kerstjens and E. Martm-Martmez, Phys. Rev. 
D 92, 064042 (2015) 

[31] E. Martm-Martmez, A. Smith, and D. Terno, (2015) 
arXiv: 1507.02688 

[32] E. Martm-Martmez, E. G. Brown, W. Donnelly, and 
A. Kempf, Phys. Rev. A 88, 052310 (2013) 

[33] E. G. Brown, W. Donnelly, A. Kempf, R. B. Mann, 

E. Martm-Martmez, and N. C. Menicucci, New J. Phys. 
16, 105020 (2014) 

[34] S. J. Summers and R. F. Werner, Phys. Lett. A 110 , 257 
(1985). 

[35] S. J. Summers and R. F. Werner, J. Math. Phys. 28, 2440 
(1987). 

[36] D. M. T. Benincasa, L. Borsten, M. Buck, and 

F. Dowker, Glass. Quant. Grav. 31, 075007 (2014) 

[37] P. Langlois, Annals Phys. 321, 2027 (2006) 

[38] J. Louko and A. Satz, Glass. Quant. Grav. 23, 6321 
(2006). 

[39] S. Schlicht, Glass. Quant. Grav. 21, 4647 (2004). 

[40] E. A. Power and T. Thirunamachandran, Phys. Rev. A 
28, 2663 (1983) 

[41] P. W. Milonni, D. F. V. James, and H. Fearn, Phys. 
Rev. A 52, 1525 (1995) 

[42] 1. Dolce, R. Passante, and F. Persico, Phys. Lett. A 355, 
152 (2006) 

]43] L. 1. Plimak and S. T. Stenholm, EPL (Europhysics Let¬ 
ters) 96, 34002 (2011), 

]44] C. M. Wilson, G. Johansson, A. Pourkabirian, 
M. Simoen, J. R. Johansson, T. Duty, F. Nori, and 
P. Delsing, Nature 479, 376 (2011). 

[45] P. D. Nation, M. P. Blencowe, A. J. Rimberg, and 
E. Buks, Phys. Rev. Lett. 103, 087004 (2009). 

[46] B. Peropadre, P. Forn-Diaz, E. Solano, and J. J. Garcia- 
Ripoll, Phys. Rev. Lett. 105, 023601 (2010) 

[47] E. Dewan and M. Beran, Am. J. Phys. 27, 517 (1959) 

[48] P. Langlois, Phys. Rev. D 70, 104008 (2004). 

[49] E. Martm-Martmez and J. Louko, Phys. Rev. Lett. 115, 
031301 (2015) 










18 


[50] J. Louko and A. Satz, Class. Quant. Grav. 25, 055012 
(2008). 

[51] R. McLenaghan, Ann. Inst. H. Poincare 20, 153 (1974). 


[52] S. Czapor and R. McLenaghan, Acta. Phys. Pol. B Proc. 
Suppl. 1 1, 55 (2008) 

[53] A. Blasco, L. J. Garay, M. Martm-Benito, and 
E. Martm-Martmez, Phys. Rev. Lett. 114, 141103 
(2015) 



